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Abstract 

Decoherence in superconducting qubits is known to arise because of a variety of environmental 
degrees of freedom. In this article, we focus on the influence of thermal fluctuations in a weakly 
damped circuit resonance coupled to the qubit. Because of the coupling, the qubit frequency 
is shifted by an amount nSvQ if the resonator contains n energy quanta. Thermal fluctuations 
induce temporal variations n[t) and thus dephasing. We give an approximate formula for the 
qubit dephasing time as a function of du^. We discuss the specific case of a flux-qubit coupled 
to the plasma mode of its DC-SQUID detector. We first derive a plasma mode-qubit interaction 
hamiltonian which, in addition to the usual Jaynes-Cummings term, has a coupling term quadratic 
in the oscillator variables coming from the flux-dependence of the SQUID Josephson inductance. 
Our model predicts that 5uq cancels in certain non-trivial bias conditions for which dephasing due 
to thermal fluctuations should be suppressed. 

PACS numbers: 
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I. INTRODUCTION 



A series of recent experiments have made it clear that it is possible to manipulate the 
quantum state of macroscopic electrical circuits based on Josephson junctions . This 

breakthrough opens the way to the realization of fundamental tests of quantum mechanics, 
up to now confined to atomic physics and quantum optics, in a solid-state physics context. 
An additional interest comes from the eventual possibility of using these circuits as building 
blocks for a quantum computer j^. In view of this latest application, it is highly desirable 
to understand better how they become dephased by environmental noise. 

To estimate the dephasing rates, the Bloch-Redfield theory assumes that the qubit is 
weakly coupled to a bath at temperature T with a memory short compared to all relevant 
timescales (white noise). In that limit, it is well known that the dephasing rate is pro- 
portional only to the low-frequency part of the environment spectral density. It becomes 
increasingly clear however that this description is inadequate in a number of cases. For exam- 
ple, the Bloch-Redfield assumptions are obviously unjustified when dephasing is due to the 
: luctuations of slow environmental degrees of freedom which typically have a 1// spectrum 
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. This is also the case when a resonance of large quality factor {Q » 1) occurs in the 
environment at a frequency comparable to the qubit, since the memory of the environment 
can not be neglected then Both processes are relevant in our experiments We 
study the quantum coherence of a circuit called the flux-qubit, measured by a DC-SQUID. 
The flux-qubit is sensitive to a number of microscopic degrees of freedom : motion of nearby 
vortices trapped in superconducting thin- films, fluctuations of the junctions critical current, 
and charge noise. In addition, it is strongly coupled to the harmonic oscillator (called HO in 
the remaining of this work) constituted by the underdamped DC-SQUID and a shunt capac- 
itor to which it is connected to improve its resolution as a detector. In recent experiments we 
observed clear signatures of the strong coupling between the two systems, manifested by the 
appearance of sideband resonances in the spectrum In the present article we investigate 
theoretically the effects of this coupling on the qubit decoherence. Thermal fluctuations of 
the photon number n stored in the oscillator shift the qubit frequency by an amount n6i>Q 
and lead to dephasing. We note that a similar effect has been recently observed in the case of 
a Cooper-pair box coupled to a waveguide resonator P, 0] . When the resonator was driven 
to perform the measurement, the qubit line was shifted and broadened due to ac-Stark shift 
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and photon shot noise. 

In the first part of this article, we propose a simple analytical formula giving the pure 
dephasing time as a function of the system parameters. In the remaining we apply this model 
to the specific case of our circuit. We start by deriving the coupling hamiltonian between 
a superconducting flux-qubit and the HO. In addition to the linear coupling term 
we find it necessary to consider the next order term which is quadratic in the oscillator 
variables. We finally investigate the dependence of the shift per photon Suq on the system 
bias parameters, namely the magnetic flux enclosed by the qubit loop and the SQUID 
bias current h- In particular, we find that it is possible to cancel the dephasing per photon 
6^0 for specific bias conditions, so that the influence of thermal fluctuations on the qubit 
should be suppressed. 



II. DERIVATION OF AN APPROXIMATE FORMULA FOR THE DEPHASING 
TIME 

Let us consider the situation where a qubit of frequency Ug is linearly coupled to an 
underdamped HO of frequency Up with a coupling strength g. The qubit is supposed to be 
an ideal undamped two-level system, whereas the HO is coupled to a bath at temperature 
T which damps its dynamics with a rate k. We are interested in the limit k << Up where 
the oscillator is underdamped. We can write the total hamiltonian as H = h[—{l/2)vqaz + 
VpO^ a + gax{a + a^)]-, where we introduced the Pauli matrices (Tx,y,z in the qubit Hilbert space 
and the usual annihilation (creation) operator a (a^) for the HO. This is the well-studied 
Jaynes-Cummings hamiltonian jlll . Il^ . Let us recall a few results useful in the following. In 
the limit where \5\ = \vq — Vp\ » g (called the dispersive regime), the energy eigenstates of 
the coupled system can be written as a function of the uncoupled energy states \i, n), where 



i = 0,1 refers to the qubit state and n to the photon state of the HO, as 
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~ \l,n) + ^ \0,n + l) 

^ _^^^^|±I|i,^) + |o,n+l) (1) 
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their energies being 



E±„ = (n+1)H±|5±/^^^^^ (2) 

One sees that because of the couphng, the energy eigenstates are shifted by a quantity 
±/i£!Ml). In the presence of n photons, the dressed qubit excited state is 1+,"^,) and the 
ground state | — , n — 1) so the qubit transition is E^^^ — E^^_^ = h{uq + 2g'^{n + l)/5). This 
means that the qubit frequency is shifted by an amount z/g„ — ^ = 'n{2g'^/S) = hSuq. 
Thus, any temporal fluctuation of the photon number will lead to dephasing |9|. Let us 
introduce the mean photon number in the HO assumed to be at thermal equilibrium n = 
1 / {explhup/ kT) — 1). The stationary photon number distribution is given by a Boltzmann 
law p{n) = (l/(fi + l))(n/(n + l))". The temporal fluctuations are characterized by the 
two-time correlation function C(r) =< n(0)n(r) >. It is possible to estimate C(r) using a 
quantum Langevin equation approach as in jl3|. We show in this way in the annex A that 



C{t) =n{n+l)exp{-K\T\) +n' (3) 

In order to quantify the effect of these fluctuations on the qubit coherence, we follow 
the analysis of Blais et al. Q]. The total phase accumulated by the qubit during a free 
evolution is (f)(t) = 2n Vq{t')dt' = 0(t) +S(f){t) where we isolated the deterministic quantity 
= 2n(i'qfi + n5uQ)t from the fluctuating 5(f){t) = 2n5uQ j^{n{t') — n)dt' . Dephasing is 
described by the quantity f^{t) =< exp{i6(f){t)) > called the dephasing factor. In the limit 
where t » the variable should have gaussian statistics so that 



U{t) = exp(- < S^ity > 12) 

= exp[-(27r5z/o)V2 /J < {n{t') - n){n{t") - n) > dt'dt"] 
Combining equations El and El we obtain that 



(4) 



f^{t) =exp[-{{27r6uo)^n{n + l)/2) [ [ exp{-K\t' - t"\) dt'dt"] (5) 

Jo Jo 



Since Jg exp{—K\t' — t"\) = {2/ K)t + {2/ K^)[exp{—Kt) — 1], we find that the long-time 
decay of the dephasing factor is given by 

1 2 
< exp(i(50(t) >= exp[--(27r5z/o)%(fi + l)-t\ (6) 
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This describes an exponential decay of time constant 



= i2n6uorn{n + 1) 

It is interesting to compare this formula with the one derived in [l^ for the case when the 
HO is driven by a coherent field of mean photon n. The photon-photon correlator is then 
C(r) = nexp {—K,\t\/2). Compared to equation |21 we notice a factor 1/2 in the exponent 
which is due to the presence of the external drive, and a replacement of the n(n + 1) by n. 
This reflects the fact that a coherent field has a poissonian distribution of photon numbers 
of variance n whereas a thermal field has a superpoissonian distribution of variance n{n + l). 

We will now show that in the limit where formula [7| applies (namely >> the same 
formula also gives the decay of the spin-echo time te- The reason is that the fluctuations of 
the photon number occur on a much shorter timescale than dephasing so that they can not 
be compensated by a refocusing pulse. More quantitatively, the phase accumulated during 
the echo sequence is 



S(j)E(t) = 2tx 5uq 
so that the fluctuations are 



t/2 nt 

{n{t') - n)dt' - / {n{t') - n)dt' 

Jt/2 



(8) 



/.t/2 i-t/2 

< 5(t)E{tf > = (27r5z/o)^[ / / {n{t') - n)(n(t") - n)dt'dt'' 

Jo Jo 

+ [ [ {n{t')-n){n{t")-n)dt'dt" 

Jt/2 Jt/2 
pt/2 pt 

-2 / {n{t')~n){n{t")-n)dt'dt"] 

Jo Jt/2 



(9) 



Obviously we need to calculate only the last term. Using the expression for the correlation 
function given earlier we find that 



t/2 pt . 

/ (n(t') - n){n{t") - n)dt'dt" = —fiin + 1)(1 - exp(-/tt/2))2 

Jt/2 



Combining with previous results we obtain 



< 6(f)E{tf >= {27i6uofn{n + 1) 



-t + —{exp{-Kt/2) - 1) - —(1 - exp{-Kt/2)f 



(10) 
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In the limit where t » k,~^, the long-time decay is still dominated by the term linear in 
t and we obtain te = t^- If on the other hand we were in the opposite limit in which the 
dephasing time is shorter than the photon correlation time, the decay of the Ramsey signal 
would be gaussian and the echo would decay with a exp(— law much slower than the 
Ramsey decay. We note that this crossover between a lorentzian and a gaussian lineshape 
when the dephasing time becomes shorter than has been observed in |j9|. 

Whereas in this reasoning we only considered the case where the qubit-HO coupling is 
linear, it can also be applied evidently for a more complex interaction hamiltonian whenever 
the qubit frequency shift is proportional to the photon number. As we will see in the next 
paragraph, this is the case for our circuit in which a flux-qubit is coupled to the plasma 
mode of its measuring DC-SQUID. 



III. QUBIT-PLASMA MODE COUPLING HAMILTONIAN 
A. Description of the system 

The fiux-qubit is a superconducting loop containing three Josephson junctions threaded 



) op com 
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by an external flux = /($o/27r) [3|, |15|, ll6[. It is coupled to a DC-SQUID detector 
shunted by an external capacitor Csh whose role is to limit phase fluctuations across the 
SQUID as well as to filter high-frequency noise from the dissipative impedance. The SQUID 
is threaded by a flux = /'('^o/2vr). The circuit diagram is shown in figure^. There, 
the flux-qubit is the loop in red containing the three junctions of phases 0j and capacitances 
Ci {i = 1,2,3). It also includes an inductance Li which models the branch inductance and 
eventually the inductance of a fourth larger junction [17^. The two inductances Ki and K2 
model the kinetic inductance of the line shared by the SQUID and the qubit. The SQUID is 
the larger loop in blue. The junction phases are called 04 and ^5 and their capacitances C4 
and C5. The critical current of the circuit junctions is written Ici (i = 1 to 5). The SQUID 
loop also contains two inductances and L2 which model its self-inductance. The SQUID 
is connected to the capacitor Cgh through superconducting lines of parasitic inductance Lg. 
The phase across the stray inductance and the SQUID is denoted (pA- The whole circuit is 
biased by a current source lb in parallel with a dissipative admittance Y{uj). Since our goal 
is primarily to determine the qubit-plasma mode coupling hamiltonian, we will neglect the 
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FIG. 1: (a) qubit biased by and SQUID biased by current 1^. (b) Simplified electrical scheme 
: the SQUID-qubit system is seen as an inductance Lj connected to the shunct capacitor Csh 
through inductance Lgh- is the flux across the two inductances Lj and Lgh in series. 

admittance Y^u). 

We start writing the total hamiltonian of the circuit shown in figure UTTXl using the circuit 
theory presented in We first choose a spanning tree containing all the capacitors as 
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shown in figure CJi. We then write the loop submatrices 

/l \ 
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We note M the mutual inductance between the qubit and SQUID loops. In the notations 



of |l8| the inductance matrices are 



( M ^ 

M La 
^ L3 
and 
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We note • the matrix elements of L^^ whose expressions can be easily computed. We 



finally obtain the total hamiltonian as 

Hs = H,,^ + (<l>o/27r)^f/(0) 

where 



(15) 



Hkin 

u{4>) 



-ELii:^cos0, + ^(Eli0.-/)^ 

+ (Ell <P^ - m2l{<pA + 05 - n - + <Pa)] + ^(04, 05, 0. 



(16) 
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The Qi {Qa) are the charges stored on the capacitors Ci (Csh) and ^(04, 05, 0^) is defined 

by 



u = +2^^(04 + 05 - fT - itM^ + - m + (04 + - (04 - <Pa?] ^^^^ 

+ i(<^5 + 0A)^ + t/60A 

Our first goal will be here to simplify this hamiltonian so that the coupling of the relevant 
degrees of freedom is made clear. We will consider here that they are only two : the qubit, in 
the two-level approximation, and the plasma mode considered, if uncoupled to the qubit, as 
a harmonic oscillator. In particular, we will neglect the SQUID junctions capacitance which 
bring additional resonances at higher frequencies, and only consider the shunt capacitance 
Csh- Our approach is justified by the fact that only the plasma mode and the qubit have 
comparable energy scales, that the plasma mode is strongly coupled to the environment 
and therefore relevant for studying dephasing and relaxation, and that it undergoes thermal 
fluctuations because of its relatively low frequency. We also observed experimental evidence 
for the qubit-plasma mode strong coupling These results are a clear indication that a 
quantum-mechanical description of the coupled "qubit-plasma mode" is indeed needed. We 
will start by doing the two-level approximation on the qubit variables. 



B. Qubit hamiltonian and two-level approximation 

The hamiltonian for the qubit alone is 

H,if, h) = H,,n - ($ o/27r)2[ V COS0, + -L( ^ 0, - ff] (18) 

We first write this hamiltonian in a two-level approximation in terms of the Pauli matrices, 
which is valid here around $o/2 because of the specific properties of the circuit eigenstates. 
We define the states |0) and |1) as the eigenstates of Hqiji, 0) = Hf \ Then when restraining 
ourselves to the 0, 1 states we have by definition Hf'' = —{hA/2)az We define Hq^\f, lb) = 
Hq{f,h)-H^q^. We have 

= (V27r)2[(l/2/n)[(E0.-/)'-(E0.-^)']] 

= ($o/27r)2[-(l/2/n)(2E0.-27r-(/-7r))(/-7r)] (19) 

= -($o/27r)2(^0,-7r)(/-7r)//n 
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FIG. 2: Qubit frequency fg as a function of the bias e for A = 5.5GHz (minimum frequency in the 
figure). The dashed hne indicates the phase-noise insensitive bias point e = where dvq/de = 

forgetting constant terms in the last equation. We now want to decompose Hg^\f,If,) 
on the |0), |1) subspace. We start by writing that Hg^\f,Ib) = (l/2)[(/ioo + ^ii)-^ + (^oo — 
hu)(yz + {hoi + hiQ)ax + (z/iqi — ihiQ)ay], where hij =< i\Hq^^\j >. For symmetry reasons 
< O|(^0j "~^)|0 >=< >= so that /loo = ^ii = 0. Indeed, the hamiltonian 

Hq'^^ is invariant under the transformation T (pi ^ — 01? 02 —02 and 03 — 2-7? — 03. 
This means that the eigenstates of Hg^^ have to also be eigenstates of T. Since = /, 
these eigenstates should be either 1 or —1 so that 0i, — 02, 27r — 03) = ±-0j(0i, 02, 03) 
and 01, — 02, 27r — 03)^ = |'0i(0i' 02; 03)P- This implies that the matrix elements 
hoo = hii = 0. 

Since we can always chose the global phases of |0) and |1) so that Hqi and hiQ are real, we 
obtain that Hj,^\f,h) = hoia^ where V = ('^'o/27r)2((0| ^ 0i|l))[-(/ - n)/lii]. We next 
define Ip = ($o/27r)(l//ii) < 0| X]0i|l e = '^Ipi^x - ^o/2)/h. In the end we can write 
the total hamiltonian as 

Hqif) = -^{Aa, + ea,) (20) 

We note that this derivation generalizes the analysis of which was made under the 
assumption that the qubit loop self-inductance is negligible. Here we retrieve the result of 
I20I 1 which showed numerically that the form of the qubit hamiltonian was little affected by 
taking into account this inductance. 
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The hamiltonian QUI yields a qubit transition frequency z/^ = -y/A^ + e^. The corresponding 
dependence is plotted in figure |21 for realistic parameters. An interesting property is that 
when the qubit is biased at e = (dashed line in figure |2I), it is insensitive to first order to 
noise in the bias variable e. 

C. Plasma mode hamiltonian 

Next we turn to the "SQUID+shunt capacitor" variables 04, 05, (pA- As already explained, 
we will here make a crude approximation and completely neglect the SQUID junctions 
capacitance. This is justified by the fact that at the bias current that we use the modes 
to which they correspond are at frequencies much higher than the qubit and plasma mode. 
Moreover, we will keep only the terms of second order in the SQUID potential, which is 
equivalent to considering the SQUID as one inductance Lj{Ib, /'). On the other hand we 
will keep in the analysis the dependence on fiux of Lj which has important effects. 

In this approximation, the only dynamical variable in the system is (pA- Its hamiltonian 
is very simply given by the hamiltonian of a harmonic oscillator of capacitance Cgh and 
inductance L = Lg + Lj{Ib, f) 

Hp = Ql/2Csh + {%/27i)\(f)A - TaY/^L (21) 

where (1)a = {'2,7i/^o)LIh is the mean value of (pA- We call a and the creation and 
annihilation operators corresponding to this harmonic oscillator : 

and is the conjugate operator. Then the hamiltonian is simply Hp = fiujQ{a^a + 1/2), 
and the phase 0^ can be written : (pA = <pA + S(po{a + a'^) (50o is the rms amplitude of the 
vacuum fiuctuations of 0a)- 

In our model where the SQUID junctions have no capacitance, for a given value of (pA 
and /' all the phases of the SQUID are well-defined functions 04,5 (0yi, /')• So the quantum 
fiuctuations of 0a translate directly into fiuctuations of 04 and 05 as follows : 

04,5 =^5 + {dMd<PA)SM<^ + + (l/2)(ci204^5/ci0^)(50o)2(a + at)2 (23) 
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where (/)4,5 = ^4,5 (0a)- We develop the functions to second order in a and a"^ for consis- 
tency. Again, the sensitivity coefficients d(j)i^5/d(j)A and (P(j)4^^^/d(j)\ depend on and /' and 
can be easily calculated. We also note that = {27i/^Q)LdIb so that 



04,5 = 0^ + ^^^Ma + a^) + i(^)'^(^'^o)^(a + (24) 

We will finally show how to evaluate the sensitivity coefficients for any SQUID parameters. 
Following 2^, Q we introduce the parameters x = (04 + 05)/2, y = (04 — 05)/2, s = 
h/{Ic4 + Ic5), b = $o/7rL5g(/c4 + /c5), Uo = $o(/c4 + /c5)/2vr, a = {Ica- Ic5)/{Ic4 + Ic5), 
and 6 = {K3 + Ki — L2 — K2)/Lsq. The stationary solutions for the SQUID phases are 
obtained by minimizing the 2-dimensional potential 

U (x, y) = Uo[—sx — cos a; cosy — a sin x sin?/ — 6sy + b{y — f /^Y] (25) 
that is to solve the equations 

dU/dx = t/o[— s + sinxcosT/ — acosxsiny] = ^^^^ 
dU/dy = Uo[cosx siny — asinx cosy — 6s + 2b{y — f /2)] = 
From this it is possible to obtain numerically the functions $4,5(0^, /') and thus all the 
sensitivity coefficients needed in the model. 



D. Qubit-plasma mode coupling hamiltonian 

The coupling hamiltonian is due to two contributions. First, the explicit coupling term 
in the equation IHEl But we also need to rewrite the plasma mode hamiltonian since the 
parameters of this hamiltonian (notably the SQUID Josephson inductance) now depends on 
the qubit state. We therefore need to reconsider the following variables in the plasma mode 
hamiltonian : 



f ^ f + (27r/<l'o)MV, 

L — ^ L + 5La,^ (27) 
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where 5L = {2TT/^o){dLj/df')MIp. The SQUID-qubit couphng term writes 

3 

V = ($o/27r)2( J] 0, - /)[/2-i'(04 + 05 - /') - /3l'(05 + <Pa)] (28) 
i=l 

Since cr^ = (X]0j ~ ^)/ < 0| we can rewrite V in the form 

y = ($o/27r)/n V,[/2-/(</.4 + 05 - /') - /3i'(05 + 0a)] (29) 
Keeping only the terms which contain exphcit couphngs we obtain that 

V = ($o/27r)/n V.[($o/27rL)(/2i'(d04/t^/6) + {I21 " l^^){dcP^I dh))6Ma + a^) 
+ {l/2){^o/2i:Lf{l^^{d?ct>,/dIl) + (/^i^ - l^r){d''ct>,/dll))5ct>l{a + at)^] 

On the other hand, the plasma mode hamiltonian now writes 



27r 



(0A - 0^)0-0: - 777? (0A ~ ^31) 



so that 



= hvoia^a + 1/2) - M|^50o(a + a^) - ' 5^l{5L/2L^){a + at) V, (32) 

Finally the total interaction hamiltonian Hj — V + Hpi_g can be written as 

Hi = h[gi{a + a^) + g2{a + at)^]^^ (33) 

The coupling constants gi and g2 could be deduced from the above expressions. Never- 
theless we propose a way to determine them experimentally. We first note that this coupling 
hamiltonian can be rewritten 



Hi = h[X,6<PA + A250i]a, (34) 

where 50^ = (pA — <f>A- This gives us a very direct way of evaluating the coupling 
constants gi and g2 : indeed varying the bias current If, through the SQUID by an amount 
SIf, is equivalent to a variation 50^ = 27rL5Ib/^o- Since we can experimentally measure the 
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dependence of the qubit frequency on the bias current e{Ib) and thus measure de/dh 
and d'^e/dll, we can obtain the couphng constants from the experiment by the following 
expressions : 

= {l/2)ide/dh)i%/2nL)5<Po 
g, = {l/A){d\/dIl){%/2'KLf6ct>l 
Finally the total qubit-plasma mode hamiltonian can be written as 



H/h = (l/2)(-Aa, - ea,) + u^ia^a + 1/2) + [giih)ia + a^) + g2{h){a + a)f]a^ (36) 

We note that a more rigorous derivation for the coupling between the plasma mode and 
the qubit after elimination of the internal dynamics of the SQUID thanks to the Feynman- 
Vernon influence functional has been performed in and gives ultimately the same form 
if we develop their interaction hamiltonian to the second order in the oscillator variab 
It is also instructive to compare this hamiltonian to the simpler situation studied in 
There a Cooper-pair box is capacitively coupled to a coplanar waveguide resonator. The 
interaction hamiltonian contains only one term, linear in the oscillator variables, and with 
a fixed coupling constant depending on the geometry of the circuit. In our case the coupling 
is mediated by the SQUID flux- dependent and current-dependent inductance ; therefore 
the coupling constants are tunable and higher-order terms are of importance. This made 
possible to induce transitions in which both the HO and the qubit state are modified 



es. 



E. Coupling constants 

We now want to give analytical formulae for the coupling constants gi and g2 in the 
simplest case where a number of assumptions are made : 1) the SQUID-qubit coupling is 
supposed to be symmetric {Ki = K2) so that the qubit bias is only coupled to the current J 
circulating in the SQUID loop e = (2/p//i)(<l>^ + MJ(4)) 2) the SQUID loop self-inductance 
and the stray inductance Lg are negligible so that the total inductance of the plasma mode 
is the SQUID Josephson inductance L = Lj{f, lb). Within these assumptions the equations 
QUI are easily solved and yield that 
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X = arcsm 



y = f'/2 



2Jccos(/'/2) 



(37) 



which imphes that the current circulating in the SQUID loop is 



J = {h-h)l2 = Ic\l 



sin(/V2) 



(38) 



,2/ccos(/'/2). 

The rms phase fluctuations of the plasma mode are 5(j)Q = (27r/$o)\/^"vW2 so that we 
obtain 
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MIp sin f'/2 



h Ale cos2 f'/2 V 2L 
MIp sin/72 hup 
h 8Ic cos2 f'/2 2L 
2 



2/CCOS/V2 



:1 -1/2 



2/ccos/V2 



2Iccosf'/2 



-3/2 



1 - 



2/ccos/V2 



-1/2- 



(39) 



Around = these expressions can be simplified by keeping only the lowest order in Ij, 



MIp sin f'/2 / hup 
h 4/c cos2 /V2 V 2L7 ' 
1 sin/V2 MIp 
16 cos"' j' 12 Lie 

We will now discuss quantitatively the behaviour of g\ and g2 for actual sample parameters 
Q : /c = 3.4M, ^ = G.5pH, Ip = 240nA, A = 5.5GHz, Up = 3.1GHz, Lj = 300pH, 
f'/2 = 1A5tc. We will restrict ourselves to a range of bias conditions relevant for our 
conditions, supposing that varies between ±300nA and that /'/2 varies by df = ±4-10^'^7r 
around 1.457r. We chose such an interval for /' because it corresponds to changing the qubit 
bias point e by ±2GIIz around 0. The constants gi and g2 are plotted in figure El as a 
function of for two different values of /' {gi is shown as a full line, g2 as a dashed line, 
and the two different values of /' are symbolized by gray for df = — 27r4 ■ 10~^ and black 
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FIG. 3: Coupling constants gi (solid line) and 52 (dashed line) as a function of the bias current, 
for two values of the reduced SQUID flux bias /' differing by df = A - 10^^ (gray and black lines). 

for df = 0). It can be seen that the couphng constants only weakly depend on the value of 

the flux in this range, so that we will neglect this dependence in the following and consider 

that gi and (72 only depend on the bias current h- Moreover we see from figure El that 

the approximations made in equation EU] are justified in this range of parameters since gi 

is closely linear in 1^ and g2 nearly constant. We also note that (71 = for Jf, = 0. This 

fact can be generalized to the case where the SQUID-qubit coupling is not symmetric and 

the junctions critical current are dissimilar : in certain conditions these asymmetries can be 

I I 

compensated for by applying a bias current for which gi{I^) = [10]. At the current 
the qubit is effectively decoupled from the measuring circuit fluctuations to first order. 

IV. ENERGY LEVELS AND DEPHASING 

In this paragraph we investigate the discuss the energy levels of the hamiltonian IHEl and 
we estimate the frequency shift per photon Suq. We will only consider the case where the 
qubit is detuned from the plasma mode Uq — Vp » gi and also |z/g — 2z/p| >> g2. 

A. Energy levels 

If the qubit and the plasma mode were uncoupled (case gi = g2 = 0) , the system energy 
eigenstates would be \i)^\n), where i refers to the qubit state and can be eother or 
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1, and n to the plasma mode occupation number. The energy levels would simply be 
E^""^ = h{iVA^ + e2+nz/o). When either gi or g2 are non zero, these eigenstates are modified, 
but for convenience we will still label them thanks to the uncoupled state from which they 
are the closest However, the energies are now modified : Ei^n = E^^^ + Sui^n- The aim 

of this paragraph is to estimate the quantity Sui^n- 

We first rewrite the hamiltonian in a more convenient way for our purpose, by introducing 
the rotated axes X and Z defined as az = (Aa^ + ea.j.)/y/ A"^ + and ax = (— ecr^ + 
Aa^) /VA2 + e2. The angle 9 is defined so that cos 9 = e/y/A^ + and sin 9 = A/VA^ + e'^. 
The system hamiltonian now writes 



H/h = -{1/2)VA'^ + e^az + iypia^a + l/2) + [gi{a + a^)+g2ia + a^f]{cos9az + sm9ax) (41) 
1. Linear term 

Let us first suppose that g2 = 0. Then the coupling is linear in the oscillator variables, 
with a longitudinal component proportional to cos 6* and a transverse component propor- 
tional to sin 6'. 

We first notice that the longitudinal component has no effect on the energy states. Indeed, 
the term az{a + a^) does not shift the energy levels to first order of the perturbation theory 
since < i, n\az{a + a'^)\i, n >= for all states. To second order, all energy levels are shifted 
by the same quantity {gicos9)'^ /up which implies that all the transition frequencies stay 
unchanged. This conclusion stays true to all orders of perturbation theory. 

On the other hand, the transverse cou plin g term gisin9ax{0' + cl^) produces the well- 
known dipersive shift in cavity QED bJ ] to second-order in perturbation theory, as 
calculated in the first section of this article. In the rotating wave approximation, we recall 
that (5z//„ = i{gi sm9y/6+{2i — l)n{gi sm9)'^/6 where 6 = a/A^ + — z/q is the qubit-plasma 
mdoe detuning. However, it is necessary here to go beyond the rotating wave approximation 
since 6 is of the same order of magnitude as the qubit and the oscillator frequency. It is 
easily seen that to second order of perturbation theory, one obtains 

Hn = 2^(^ismg)^ J^/_^, + {2^ - l)2n{g, sin 9 f ^^^^ (42) 
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The first term of tliis equation describes tfie Lamb sliift, wliicli simply renormalizes tfie 
bare qubit frequency and fias no influence on dephasing. We will thus neglect it in the 
following. The frequency shift per photon is 



6u^ = 4{gi sin 9) 



2 VA^ + e^ 
A2 + e2 - 1/2 



(43) 



From the previous expression it is clear that the sign of Suq is fully determined by the 
sign of 6. In particular, in our experiments where Uq > Up, 5z/q > 0. 

2. Quadratic term 

Next, we consider the case gi = but g2 > (which is the case notably at the decoupled 
current h = The quadratic coupling term produces effects which are sensibly 

different from the cavity QED case. Indeed, it generates a frequency shift to first order in 
perturbation theory via the term 2g2 cosOazo^ a. Considering that the g2 coupling is already 
second order, we only keep the first order of perturbation theory. We therefore obtain that 
^^i,n — — (2^ — 1)2(72^ cos 6' so that the shift per photon is 



Contrary to the shift produced by the linear coupling term, the sign of this frequency 
shift now depends on e. Since g2 is negative (see figure Ej), 5^1 actually has the same sign 
as e. We also note that the quadratic term has no effect on the qubit when e = 0, since at 
that point the average flux generated by both qubit states |0) and |1) averages out to zero 
so that the SQUID Josephson inductance is unchanged. 

3. Total frequency shift and dependence on the bias parameters 

The total frequency shift per photon is the sum of the two contributions identified above 




4(72 cos 9 



(44) 



5iJo{e,h) = Suq + 5ijI 



4 {gismOf 



v/A2 + e2 



- - 5f2 COS 9 



(45) 



A2 + e2 - z/2 
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FIG. 4: Frequency shift per photon 5i'q as a function of and e. The white regions correspond 
to —15MHz and the black to +35MHz. The sohd hne €m{Ib) indicates the bias conditions for 
which di^Q = 0. The dashed hne indicates the phase noise insensitive point e = ; the dotted hne 
indicates the decouphng current If, = I^. 

Because of the different dependence on e of the two contributions discussed above, we 
expect a cancellation of the AC-Zeeman term (due to gi) by the quadratic term (due to (72) 
for some bias parameters corresponding to a negative value of e. This is shown in figure E] 
where we plotted Suole, h) as calculated with the formula above for the sample parameters 
considered in the previous paragraph. The curved full line corresponds to the points em{h) 
for which Suq = 0. For these bias conditions, it is expected that the qubit is insensitive 
to the thermal fluctuations of the plasma mode (see formula [7j). Therefore we predict an 
increase of the dephasing time whenever e = em{h)- 

We stress that these biasing conditions are non-trivial in the sense that they do not satisfy 
an obvious symmetry in the circuit. This point is emphasized in figure |3] where we plotted 
as a dashed line the bias conditions e = for which the qubit is insensitive to phase noise 
(due to flux or bias current noise) ; and as a dotted line the decoupling current conditions 
lb = II for which the qubit is effectively decoupled from its measuring circuit. The em{h) 
line shares only one point with these two curves : the point (/^,e) which is optimal with 
respect to flux, bias current, and photon noise. For the rest, the three lines are obviously 
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distinct. This makes it possible to experimentally discriminate between the various noise 
sources limiting the qubit coherence by studying the dependence of on bias parameters. 



V. CONCLUSION 



Superconducting qubits are often measured by circuits behaving as underdamped os- 
cillators to prevent energy relaxation of the qubit. If these oscillators have a frequency 
comparable to kT, their photon number undergoes thermal fluctuations. This induces fre- 
quency dispersive frequency shifts of the qubit frequency uSuq and leads to dephasing. In 
this article we derive a simple fomrula to account for this process. We apply our model 
to the specific case of a fiux-qubit coupled to the plasma mode of its DC-SQUID. Because 
of the SQUID internal degrees of freedom (circulating current), the interaction hamiltonian 
contains two terms, one linear in the oscillator variables which describes an effective induc- 
tive coupling between the two circuit, but also a quadratic term due to the fiux-dependence 
of the SQUID Josephson inductance. Moreover, the coupling constants can be tuned over 
a wide range by changing the SQUID bias current. We study the qubit frequency shift per 
photon Suo and find that Suq = for non-trivial biasing conditions. When they are fulfilled, 
the effect of thermal fluctuations on the qubit should be suppressed. 



VI. ANNEX A 



Here we show how we evaluate the correlation function C{t) =< {0)a{0)a^ {t)a{t) >. In 
order to do so, we follow ^|. We model the damping of the HO by a linear coupling to a 
bath of harmonic oscillators 



H — Hho + Hi^ath + Hint , 

(46) 

= + I) + nuj{A\Lo)A{Lo) + i) + h\g{Lo)a^A{Lo) + g*{u)A\Lo)a] 

Under the assumption that the bath has a short memory, it can be shown that the 
evolution of the HO variables in the Heisenberg representation is given by 



a= {-i2TiPp- K/2)a{t) - Fit) (47) 
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where F(t) = g{u)A{u!,0) exp{—iut) is a quantum-mechanical operator describing the 
random force acting on the HO. The damping rate of the field, which we write k/2 since in 
our notations n is the damping rate of the intensity stored in the HO, can be related to the 
bath parameters. Introducing the density of states //(tu) we have k/2 = iTri{2TTi'p)\g{27Ti'p)\'^. 
Integrating equation HTJ we obtain 



a{t) = a(0) exp (-i27rz/p - K/2)t - exp (-?27rz/p - K/2)t / F{t') exp {i2TTUp + K/2)t'dt' (48) 



C{t) = < at(0)a(0)at(t)a(t) > 

= < n(0)^ > exp(-Kt) 

- exp {-Kt) /g < at(0)a(0)Ft(t')a(0) > exp (-i27rz/p + (/t/2))t' 

+ exp {-Kt) < at(0)a(0)at(0)F(t') > exp {i27TUp + (/«/2))t' 

+ exp i-Kt) /g /J < at(0)a(0)Ft(t')F(t") > exp [-i27iUp{t' - t")] exp [(/€/2)(t' + t")]dt'dt") 



In this equation, the second and third term vanish. Indeed, 

< at(0)a(0)Ft(t')a(0) > = ^ di^) exp"'"* < a^{0)a{0)a^{0)A{uj, 0) > 

(50) 

= ^Ec^fHexp"*^* < at(0)a(0)at(0) >< A{uj,0) > 
since it is assumed that at time t = the bath and the HO are uncorrelated. For a bath 
at thermal equilibrium, < A{uj,0) >= so that < a'^{0)a{0)F'^{t')a{0) >= 0. The same 
reasoning holds of course to show that < aUO)a(O)aUO) F (f) >= as well. Using the fact 
that < F\t)F{t") >= KN{2TxVp)5{t' - t") [13], where N{uj) =< A^{uj,0)A{uj,0) >, we can 
calculate the last term : 

exp i-Kt) /J /q < at(0)a(0)Ft(t')F(t") > exp [-i27TUp{f - t")] exp [(/t/2)(t' + t")]dt'dt" 
= < n{0) > kN{2tiVp) exp {-nt) j'l j'l 6{t' - t") exp [-i2'KVp{t' - t")] exp [^(t' + t") /2]dt'dt' 
= < n(0) > KN{27iUp) exp (-Kt) /J exp {Kt')dt' = N{2TCUp) < n{0) > (1 - exp (-/tt)) 



Since we assume that the HO and the bath are permanently in thermal equilibrium, 
N{2TiVp) =< n{0) >= n, whereas < n^(0) > — < n{0) >^= n{n + 1) (non-poissonian 
photon statistics of a thermal field) . Therefore we obtain 




This allows us to calculate C{t) =< {0)a{0)a'^ {t)a{t) > 



(49) 



(51) 
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C(t) = [< 71^(0) > - < n(0) >^] exp {-kt)+ < n(0) >^ 
— n{n + 1) exp {—kt) + 



(52) 
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